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Abstract-This paper presents an efficient method for the de- 
sign of two-channel quadrature mirror filter (QMF) bank by 
polyphase representation of low-pass prototype filter and power 
optimization method for minimization of an objective func- 
tion. The objective function is formulated as linear combina- 
tion of pass-band error and stop-band residual energy of the 
low-pass analysis filter of the filter bank. Our simulation 
results show that the use of polyphase components in the 
design of low-pass prototype filter reduces the computational 
complexity and the number of iterations. Two design examples 
are included to show the validity of the proposed method and 
its superiority over existing methods. 

Keywords — Filter banks, Polyphase component, Power 
optimization method, QMF bank, Eigenvector. 

I. Introduction 

A two-channel QMF bank is extensively used in many 
signal processing fields. Initially, it was used for aliasing 
cancellation in sub-band coding of the speech and image 
signals [1-4]. Now a days the applications of QMF can be 
found in design of wavelet bases [5,6], speech and image 
compression [7,8], trans-multiplexers [9], antenna systems 
[10], digital audio industry [11], and biomedical signal 
processing [12] due to advancement in QMF banks. Because 
of such wide applications, many researchers giving a lot of 
attention in efficient design of such filter bank. 

A typical two-channel QMF bank is shown in Fig. 1, the 
reconstructed signal y(n) is not perfect replica of the input 
signal x(n) due to three types of errors: aliasing distortion 
(ALD), phase distortion (PHD), and amplitude distortion 
(AMD). ALD can be cancelled completely by selecting the 
synthesis filters cleverly in terms of the analysis filters, 
whereas PHD eliminated using the linear phase FIR filters 
[1,17]. AMD is then minimized by optimizing the filter tap 
weights of the low-pass analysis filter using computer aided 
techniques. Several techniques for the design of QMF bank 
have been reported [8, 13-24]. 

The design problem for the QMF bank is to find the filters 
of the analysis/synthesis sections such that the reconstructed 
signal y(n) approximates the original signal x(n) in some 
optimal sense. 

The z-transform of the output signal y(n) of the two-chan- 
nel QMF bank, is given by [ 1 3-1 5, 1 9] . 
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Fig. 1 Two-channel quadrature mirror filter bank 

%[tf„(-z>fb 00 + H t (-z)F t (z)l X(-z) (1) 

aliasing can be completely removed by defining the synthesis 
filters as given below 

F Q0 = #i£-zj and F L (z) = -H D (-z) (2) 

The relationship H;{z) = i?;(-z) between the mirror image 
filters and by using "(2)," the expression for the alias free 
reconstructed signal can be written as: 

Viz) = lH (z} Rii-z) - H L G:)fi D (-z>]^Cz} 

= ¥A.Rliz')-Rt{-zy\Xiz-) ( 3 ) 
or ¥(z} = T(z-)X(z) (4) 

If the analysis filter H (z) is selected to be a linear phase FIR 
then phase response of the transfer function , T(z) , of the 
QMF bank also become linear and phase distortion of the 
QMF bank is eliminated. The corresponding frequency re- 
sponse can be written as 

ffnC*^*) = B-^ N - L ^ l H R (ca} (5) 
where H R (a>) is the amplitude function. Hence, by substituting 
"(5)," the overall transfer function of the QMF bank can be 
written as 

Ite-^^mJ**?? - (-ir- l \n y&-»)\*] 

(6) 

If the length of filter, N, is odd, above equation gives T(e> m ) = 
at co = n/2, implying severe amplitude distortion. Therefore, 
N must be chosen to be even to avoid this distortion. By 
"(6)," the condition for perfect reconstruction is given as 

IH^W + IH^"-")? =c (7) 
Now by "(3)," the design problem of filter bank reduces to 
the determination of the filter tap coefficients of the linear 

-ACEEE 



Short Paper 



ACEEE Int. J. on Signal & Image Processing, Vol. 4, No. 1, Jan 2013 



phase FIR low-pass analysis filter H g (z) only, subject to the 
perfect reconstruction condition of "(7),". 

This paper proposes an improved method for the design 
of a two-channel QMF bank in which polyphase components 
are used for implementation of filter bank. The polyphase 
component based design of QMF bank is efficient in the 
sense of speed constraints imposed on the processors. The 
objective function, which is a linear combination of the pass- 
band error and stop-band residual energy of the low-pass 
analysis filter HJz) is minimized using power optimization 
method. The results are compared with the other existing 
methods to demonstrate the performance of the proposed 
method. 

The organization of the paper is as follows. Section 2 
gives the brief description of the implementation of QMF 
bank through polyphase component. Section 3 describes the 
design problem formulation. Section 4 describes the proposed 
optimization method and design algorithm. Section 5 examines 
the performance of the filter bank designed and comparison 
with other existing methods. 

II. QMF Bank Using Polyphase Component 

The QMF bank can be realized efficiently by using 
polyphase structure [17], which enables us to rearrange the 
computation of the filtering operation. For decimation filter, if 
direct form implementation is used then only the even num- 
bered O/P samples are computed and this computation re- 
quires (N+l) multiplication per unit time (MPUs) and N addi- 
tion per unit time (APUs). However during the computation 
of odd numbered O/P samples, the structure is merely rest- 
ing. If we use polyphase implementation then the computa- 
tion of O/P samples requires (N+ 1)/2 MPUs and N/2 APUs 
[1]. The multipliers and adders in each of the filter E Q (z) and 
E (z) now have two units of time for doing their works and 
they are continually operative and no resting time. Thus 
polyphase representation of decimation filter, as shown in 
Fig. 2, reduces the computational complexity of the multipli- 
ers and adders in filter bank. 
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Fig. 2 Polyphase implementation of decimation filter 

The prototype filter HIz) can be written in its polyphase 
form as 



and 



by putting the value of £,(z) in to "(8)," we get 
H„Ge> = E (z 2 } + z- N E (z- 2 ) 



or 



(8) 
(9) 

(10) 
(11) 



where 

= 3£;r i3 ''% We-*™ (12) 

and = SfiT^Si (n-) e V™ (13) 

so, the equation of prototype filter can be written as: 

I? =D e D + ^ (14) 

The amplitude response of prototype filter can be given as: 

or 

d5) 

2 j'q>JV j- . r ■. 

1^(^)1 =l£? =D Ms" ! cos(«(2n-y)}| 2 



T T 
= p cc p 

where vectors p and c are: 
P = lPo Pl Pi 



■PmY 



(16) 



(17) 



c = [cos bt (— ~) cosiri {2 — -^j ,„ .- cos isi (lM — ~j] T 

(18) 

The vectors p and c will be used for the designing of objective 
function in the next section. 

III. Formulation of Design Problem 

The objective function <p to be minimized for the design 
of QMF bank using polyphase components can be chosen 
as a linear combination of two functions: 



(19) 



The constant a is real and in the range 5 a 5 1 , E and 
E are the measure of stop band residual energy and pass 
band error respectively. We propose to minimize the above 
objective function, by optimizing the coefficients of low pass 
filter H (z). The calculation of stop band error and pass band 
error is given below. 

A. Pass band Error (E ) : 

The amplitude response at zero frequency is given by 

H K (0j = p T l (20) 

where 1 is the vector of all l's. By taking this as a reference 
then the pass band deviation at any frequency can be given 
as: 

p T l - p T c = p T (l - c) (21) 
where vector c is defined in "(18)," so the pass band error 
can be written as: 



e f = P T QP 
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where 

Q= i/^Cl-cXl-tf^ (23) 

B. Stop band Error (E ): 

The stop band error E is given as: 

£,= ij-;jif(Oj=^= p T Sp (24) 

where 

it J -m i 

vector c is defined in "( 1 8)," and by using the value of E and 
£ the objective function can be written as : 

<p = p T R p (26) 

where 

R = o-S + Cl - oQQ (27) 

Here matrix R is a real, symmetric and positive definite matrix. 
The unit norm vector p which minimizes the objective 
function <p is the eigenvector corresponding to the minimum 
eigenvalue, I of matrix R which is calculated using 'Power 

G nun, ° 

method' described in next section. 

IV. Power Optimization Method and Design Algorithm 

The power method [1] is a simple iterative method which 
is designed to compute the dominant eigenvalue and 
corresponding eigenvector of the matrix. The other standard 
methods for finding the eigenvalues of a matrix R are 
impracticable when R is large. Even evaluating the determinant 
of an n x n matrix is a huge task when n is large. The principle 
of power method is that the vectors multiplied by matrix R, 
stretched in the direction of eigenvalue. In this method, we 
start with an initial vector h Q chosen arbitrarily and generate 
a sequence of approximations to get optimized filter 
coefficients. Here the iterative steps are performed a total of 
n times unless the vector h Q is orthogonal to every vector in 
the eigenspace. The matrix R contains complete set of 
eigenvectors with corresponding eigenvalues. This 
eigenvector corresponding to minimum eigenvalue contains 
the final optimized filter coefficients. 

Design Algorithm: 

Step 1: Guess an eigenvector h o and normalize it as: 

y° = h %h D ii 

Step2: Assume the length of the filter ( N). 
Step 3: Calculate inverse of matrix R and put it equal to 
R 1 . 

Step 4: Multiply y by R t to get a new vector. Normalize 

the result and call it as y l _ y . R r 
Step 5: Repeat the process n times so that 

%3*-+ = = *$a 

Hence generated vector y n is the eigenvector of matrix 'R' 
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which minimizes the objective function and 1/1 is the 
minimum eigenvalue (A ) of matrix 'R' . 

V. Results and Discussion 

A MATLAB program has been written to implement the 
design procedure described in the previous section. Use of 
polyphase components with aliasing free condition makes 
the system more efficient by reducing the computational time 
and number of iterations (NOI). The comparison of the 
proposed method for the design of QMF bank with other 
existing method is given in Table 1 . Two design examples are 
presented to illustrate the effectiveness of the proposed 
method. 

The performance is evaluated in terms of pass band error 
(E), stop band error (£), NOI, computational time (CPU time), 
stop-band first lobe attenuation (A ; ), stop-band edge 
attenuation (A) = " 20 log 10 (H ((O s )) as given in [ 19] and peak 
reconstruction error (PRE) in dB. The filter coefficients initially 

taken are h Q (n) = [0, 0, 0, 0, 0.707]. 

Example 1: For filter length (N) = 24, to = O.671, co = 0.471 and 
oc = 0.52. 

The optimized 12 filter coefficients of low pass analysis filter: 
h(n) = [-0.0025 0.0173 -0.0001-0.0361 0.0076 0.0628 

-0.0287-0.1075 0.0771 0.2006 -0.2556 -0.9330] 
The corresponding magnitude plots of analysis filters 
H g (z) and H^z) are depicted in Fig. 3a. The attenuation 
characteristics of low-pass filter H Q (z) is plotted in Fig. 3b. 
Fig. 3c and 3d represent the magnitude of distortion function 
T(z) and reconstruction error of QMF bank (in dB) 
respectively. The significant parameters obtained are E = 
2. 1975x10 s , E = 2.924x 10 \ A = 26.84 dB, A, = 36 dB, NOI=02, 
PRE = 0.2249 dB, and CPU-time = 0.00295sec (on core 2 duo 
processor 2. 1 GHz, 1 GB RAM). 
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Fig. 3a Amplitude response of analysis filters for N =24 
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Fig. 3b Attenuation characteristics of low-pass analysis filter 
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Fig. 3c Amplitude of the distortion function 
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Fig. 3d Reconstruction error in dB 

Example 2: For filter length (N) = 32, eo = 0.6n, co = OAn and 
a = 0.72. 

The optimized 16 filter coefficients of low pass analysis filter: 

h(n) = [-0.0005 0.0062 0.0011 -0.0140 -0.0004 0.0255 

-0.0041 -0.0430 0.0141 0.0682 -0.0360 -0.1096 
0.0850 0.1998 -0.2620 -0.9291] 
The corresponding magnitude plots of analysis filters 
H g (z) and H^z) are shown in Fig. 4a. The attenuation 
characteristics of low pass filter H Q (z) is plotted in Fig. 4b. 
Fig. 4c and 4d represent the magnitude of distortion function 
T(z) and reconstruction error of QMF bank (in dB) 
respectively. The significant parameters obtained are E = 
4.227xl0 7 ,£ = 1.98x10 s , A= 37 dB, 
A = 45.86 dB, NOI = 03,' PRE = 0.1515 dB, and CPU- 
time = 0.00663sec. 




Fig. 4a Amplitude response of analysis filters for N=32 

Table I summarizes that the proposed method gives 
improved performance than all other existing methods in terms 
of number of iteration (NOI) required and the computational 
time of the processor. Therefore the proposed method is 
computationally efficient. 
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Fig. 4c Amplitude of the distortion function (JV=32) 
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Fig. 4d Reconstruction error in dB (N=32) 

Table I. Performance Comparison of Proposed Method with Other 
Existing Methods for N =32 
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Conclusion 

In this paper, a new technique for designing of quadrature 
mirror filter (QMF) bank using polyphase component has 
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been presented. The design problem minimizes a weighted 
sum of pass band error and stop band residual energy of the 
low-pass analysis filter. The comparison of the simulation 
results indicates that the proposed method requires less 
computational time ( few msec, on core 2 duo processor 2. 1 
GHz, 1 GB RAM ). The Polyphase decomposition results 
dramatic computational efficiency and reduces the number 
of iterations. This approach can be extended to the case of 
two-dimensional filters and for filters with larger taps. 

References 

[1] Vaidyanathan, P.P. (1993). Multirate systems and filter banks. 

Prentice Hall, Englewood Cliffs, NJ. 
[2] Smith, M. J. T., Eddins, S. L. (1990). Analysis / synthesis 

techniques for sub-band image coding. IEEE Trans. Acoust. 

Speech Signal Process., ASSP-38(8), 1446-1456. 
[3] Vetterli, M. (1984). Multidimensional sub-band coding: Some 

theory and algorithms. Signal Process. 6, 97-112. 
[4] Liu, Q. G., Champagne, B., Ho, D. K. C. (2000). Simple design 

of oversampled uniform DFT filter banks with application to 

sub-band acoustic echo-cancellation. Signal Process. 80(5), 

831-847. 

[5] Chan, S. C, Pun, C. K. S., Ho, K. L. (2004). New design and 
realization techniques for a class of perfect reconstruction 
two-channel FIR ulter banks and wavelet bases. IEEE Trans. 
Signal Process. 52 (7), 2135-2141. 

[6] Sablatash, M. (2008). Designs and architectures of filter bank 
trees for spectrally efficient multi-user communications: review, 
modifications and extensions of wavelet packet filter bank 
trees. Signal, Image and Video Processing (SIViP), Vol. 5, No. 
1, 09-37. 

[7] Xia, T., Jiang, Q. (1995). Optimal multiulter banks: Design 
related symmetric extension transform and application to image 
compression. IEEE Trans. Signal Process., 47 (7), 1878-1889. 

[8] Kumar, A., Singh, G. K., Anand, R. S. (2010). An improved 
method for the designing quadrature mirror filter banks via 
unconstrained optimization. /. Math. Model. Algorithm 9(1), 
99-111. 

[9] Bellanger, M. G, Daguet, J. L. (1974). TDM-FDM trans- 
multiplexer: Digital polyphase andFFT. IEEE Trans. Commun. 
22 (9), 1199-1204. 

[10] Chandran, S. (2003). A novel scheme for a sub-band adaptive 
beam forming array implementation using quadrature mirror 
filter banks. Electron. Lett., 39(12), 891-892. 



[11] Painter, T., Spanias, A. (2000). Perceptual coding of digital 

audio. Proc. IEEE, 88(4), 451-513. 
[12] Afonso, V. X., Tompkins, W. J., Nguyen, T. Q., Luo, S. (1999). 

ECG beat detecting using filter banks. IEEE Trans. Biomed. 

Engg., 46(2), 192-202. 
[13] Chen, C.K., Lee, J.H. (1992). Design of quadrature mirror 

filters with linear phase in the frequency domain. IEEE 

Trans.Circuits Syst. 39 (9), 593-605. 
[14] Jain, V. K., Crochiere, R. E. (1984). Quadrature mirror filter 

design in time domain. IEEE Trans. Acoust. Speech 

Signal Process. ASSP-32 (4), 353-361. 
[15] Xu, H., Lu, W. S., Antoniou, A. (1998). An improved method 

for the design of FIR quadrature mirror image filter banks. 

IEEE Trans. Signal Process.46 (6), 1275-128 1 . 
[16] Nayebi, K., Barnwell in, T. P., Smith, M. J. T. (1992). Time 

domain filter analysis: A new design theory. IEEE Trans. 

Signal Process. 40 (6), 1412-1428. 
[17] Vaidyanathan, P. P. (1 990). Multirate digital filters, filter banks, 

polyphase networks and applications: A tutorial. Proc. IEEE 

78 (1), 56-93. 

[18] Lu, W. S., Xu, H., Antoniou, A. (1998). Anew method for the 
design of FIR quadrature mirror-image filter banks. IEEE Trans. 
Circuits Syst. II: Analog Digital Signal Process. 45(7), 922- 
927. 

[19] Sahu, O. P., Soni, M. K., Talwar, I. M. (2006). Marquardt 
optimization method to design two channel quadrature mirror 
filter banks. Digital Signal Process. 16(6), 870-879. 

[20] Kumar, A., Singh, G. K, Anand, R. S. (2011). An improved 
method for the design of quadrature mirror filter bank using 
the Levenberg-Marquardt Optimization. Signal, Image and 
Video Processing ( SIViP ). 

[21] Upendar, J., Gupta, C. P., Singh, G K. (2010). Designing of 
two channel quadrature mirror filter bank using Partical Swarm 
Optimization. Digital Signal Processing 20(10), 304-313. 

[22] Kumar, A, Singh, G. K, Anand, R. S. (2009). A closed form 
design method for the two-channel quadrature mirror filter 
banks. Signal, Image and Video Processing (SIViP), Vol. 5, 
No. 1, 121-131. 

[23] Gupta, A., Agrawal, S. K. : Designing of Two Channel 
Polyphase Quadrature Mirror Filter Bank using Power 
Optimization Method. International Conference on Computer 
& Communication Technology (ICCCT), (2011) 

[24] Sahu, O. P., Soni, M. K, Talwar, 1. M. (2006). Designing 
quadrature mirror filter banks using steepest descent method. 
Jorunal of Circuits Systems and Computers, Vol. 15, No. 2, 
29-42. 



©2013 ACEEE 
DOI:01.IJSIP.4.1.13 



28 



-iV ACEEE 



